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ABSTRACT 

I consider the recent discovery of a soft X-ray source inside the error box of the gamma ray 
burst GB 960720 by the SAX, ASCA and ROSAT satellites, in terms of the fireball model. I 
show that the ejecta shell, which, after causing the burst is cold and dense, but still relativistic, 
keeps plowing through the interstellar medium, heating up the just-shocked matter which 
then emits X-rays. I compute the radiation emitted by this matter. I show that, up to about 
two months after the burst, in the cosmological scenario a soft X-ray (0.1 — 10 keV) flux of 
at least ~ 10“^^ erg cm~^, well within current observational capabilities, is generated, 
explaining the observations of the three satellites. Instead, in the Galactic Halo scenario a flux 
3 orders of magnitude lower is expected. Detection of this non-thermal, declining flux in a 
statistically significant number of objects would simultaneously establish the fireball model and 
the cosmological nature of gamma ray bursts. 


Subject headings: gamma rays: bursts - radiation mechanisms: nonthermal - X rays: general - 
hydrodynamics - relativity - shock waves 


1. Introduction 

Our current theoretical understanding of gamma ray bursts (GRBs) is mostly based upon the fireball 
model (Meszaros, Laguna, Rees 1993). However, though this model has won great critical acclaim because 
of its ability to explain two otherwise mysterious features (the bursts’ time duration and the nonthermal 
spectra), it still has made no testable predictions that would allow gauging it against observations. This is 
mostly due to our ignorance of the physics of electron acceleration at relativistic shocks. To get a feeling of 
how serious this problem is, notice that the very same, simplified analysis of the radiation emitted behind 
the shock was carried out in two papers treating exactly the same physical problem, but separated by 
twenty years, Blandford and McKee (1977, dealing with AGNs) and Sari, Narayan, Piran (1996, dealing 
with GRBs), and that furthermore the first one was published before particle acceleration at shocks was 
even discovered (Bell 1978). 

It is the aim of this paper to derive a prediction from the fireball model, by (nearly completely) 
circumventing the problem of electron acceleration at relativistic shocks. Surprisingly, the prediction 
is different for cosmological and Galactic Halo scenarios. The opportunity to do this is offered by the 
detection, by the X-ray satellite SAX (Piro et ai, 1995), of the gamma ray burst GB960720, both in the 
hard X-ray/soft y-ray band, where its results are confirmed by simultaneous observations by BATSE in 
a nearly identical band, and in the soft X-ray band (Piro et ai, 1996a). Subsequent observations of GB 
960720 (Piro et ai, 1996b, Murakami et ai, 1996, Greiner et al.^ 1996) made ~ 45 days after the burst. 
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have shown that a weak source, which is not an AGN, is present inside the WFCs’ error box, 5 arcmin, 
with a flux « 2 X 10“^^ erg s~^ cm~^ in the band 0.1 — 10 keV. 

It is thus interesting to speculate about what a systematic search for GBRs’ afterglow ought to yield, 
in the fireball model. This afterglow has a simple interpretation: it is the cooling of matter swept up 
by the (still relativistic) shell of ejecta, plowing through the interstellar medium. In the next Section, 1 
compute the expected soft X-ray fluxes from GRBs some time after the burst, for the fireball model, both 
in the cosmological and in the Galactic Halo scenarios. In the last Section, I discuss how the theoretical 
computations relate to the SAX, ASCA and ROSAT observations, and why the soft X-ray band is ideal for 
carrying out a statistical search for a soft X-ray afterglow. 


2. Predictions for the cosmological and Galactic Halo scenarios 


In the fireball model (Meszaros, Laguna, Rees, 1993), an initial energy E is released, together with a 
contaminating mass of baryons M^j = E/rjc^, with rj ^ 10^, as required by observations. In the cosmological 
scenario, E = erg, and M^j ^ 10^® g. For reasons to be explained later, I shall concentrate on 

long 1 s) bursts, for which 7 100 (Sari and Piran 1995). In this case, after a phase of free expansion, 

when the expanding ejecta have swept up an ISM mass « Mej/g Mgj, a shock with the ISM forms, at a 
distance from the site of energy injection 

Rsh = 5 X eK^ (^) 

where n = 1 ni cm~^ is the ISM density. At this point, a reverse shock propagates backward toward the 
still freely expanding ejecta, converting their directed kinetic energy into internal energy, whose prompt 
release (on a timescale of seconds) leads to the GRB. 

Gonventional analysis halts here, but the shell still has a Lorenz factor g/2 (Sari and Piran 1995), so 
that it will continue to expand relativistically into the interstellar medium. Post-shock material has cooling 
times of order 10® s in the shock frame (Meszaros, Rees, Laguna 1993, Sari, Narayan, Piran 1996), 
depending on various details, but always negligible with respect to hydrodynamical evolution times of the 
post-burst shell which, as will be shown later, are of order « 1 month. Thus, the shell can be regarded as 
cold and dense; its evolution is given by Blandford and McKee (1976) as 

'y ~ ^ ~ i\,p i)+2 ~r (2) 

MT r]-2 

ej ' 

where I used the fact that the shell has initial Lorenz factor g/2. This equation is exact for any shock 
(and shell) speed, sub or super relativistic. Here M is the present mass in the shell, M = Mej + nrripV 
where V is the volume swept up. From the above equation, we see that the shell is relativistic (7 ^ 2) until 
M Ki y/3Mej, i.e., till it reaches a radius 

Rrei ~ rj^^^Rsh = 2.4 X (^) 


As long as the expansion is relativistic, eq. g can be rewritten, defining an adimensional time since the 
burst X = ct/Rsh, as 


7 = 1 + 




( 4 ) 





- 3 - 


Since the cooling time is much shorter than the shock evolutionary time, the total energy radiated per unit 
time, a relativistic invariant, also follows from hydrodynamical arguments (Blandford and McKee 1976): 

E = A'KR^v"i{'^ — l)nmpC^ (5) 


where R, 7 and v are the shock position, Lorenz factor and speed, respectively. It is convenient to rewrite 
Eq. H as 

( 1 on \ ^ 

— j fix) erg , (6) 

where 

fix) = (1 + x)^(7 - 1)^7^ - 1 . (7) 

Because of the superluminal expansion effect, the time as measured by an observer on Earth, ^ 0 , is given 
by dt^ = (1 — V/c)dt] I define an adimensional Earth time as 


dy = 


c dt^ 

Rsh 


(1 _ 


( 8 ) 


Eqs. ^ and together define parametrically the dependence of the luminosity in terms of Earth time; this 
is plotted in Fig. 1, for the range of time for which the shell is relativistic. From this (and Eq. ^ we 
see that the shell is relativistic for « 2 month after the burst. From Eq. ||, it can be seen that the total 
distance (and thus the time) before slowdown to sub-relativistic speed is only mildly dependent upon y 
(like vf-^^), so that the total lapse of time before afterglow turnoff is reasonably well determined and, as 
pointed out to me by the referee, roughly the same for both short and long bursts, which are thought to 
differ mostly because of their different ry-values. Also, from Fig. 1 it can be seen that /(y) is a curve with 
some curvature, so that no power-law dependence of E on time t is meaningful. Lastly, notice that the 
exact shape of Fig. 1 is somewhat sensitive to rj, though not its normalization after « 1 month. 


In the cosmological scenario, for a typical source distance D = 1 Gpc, I obtain for the total flux 
radiated by the shell in the relativistic snowplow phase, from Eq. g 


F = 10 erg s' 


-1 -2 l/3c,2/3 

cm n.^ E^[ 


100'N 
h ) 


4/3 


I Gpc 
D 


fiv) 


( 9 ) 


from which we see that, in the whole first month after the burst, the total flux level always exceeds 
« 3 X 10“^^ erg cm^. 

Before estimating which fraction of this flux ends up in the soft X-ray regime, I derive the total flux 
in the Galactic Halo scenario. Since in this case sources are « 10^ times nearer, the total energy released 
is 10® times lower, E = 10'^® erg. Since approximately the same value of 77 = 10^ is required both by burst 
duration and by the spectra, the contaminating baryon mass is now Mf.j = E/ijc^ = 10^*^ g (Begelman, 
Meszaros, and Rees, 1993). In this scenario, a reasonable ISM density is ni « 10“®. 

To estimate the total luminosity after about one month, I find that the total distance covered by the 
shell, if it were still relativistic, would be i? = lO®"^ cm, containing a mass Mgw = 2ni x 10^^ g. From Eq. 
Eq. H we see that the shell cannot be relativistic, because 1. 

Rather than integrating Eq. || exactly, the following simple argument will be used. The shell ceases 
to be relativistic for R^ei ~ V^^^Rsh as before, and, from Eq. |^, I find Rrei ~ I-l x 10^® cm. This occurs 
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3.5 X 10^ s after the burst. From that moment on, the expansion is subrelativistic and the distinction 
between t and % ceases to be important. In this limit, Eq. || becomes (Blandford and McKee 1976) 


Mv = 2MejC 



2MejC , 


( 10 ) 


which is the well-known nonrelativistic snowplow model, for which the shock position and speed scale as 
R cx V (X t Using as initial values the previously determined position and time at which the 
shell becomes subrelativistic, I find that, a month after the burst, the shock speed is u « 10® cm s“^, and 
R = 3.5 X 10^^ cm. With these values the total flux emitted can be computed from the subrelativistic limit 
of Eq. H, U = AirR'^v^nmp. Using a typical source distance of 100 kpc, I find 


F = 1. X 10 erg s ^ cm ® 


_2 /3 X 10® AOO kpc 


\ D 


10 


-3 


-3 \ 1/4 
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for the Galactic Halo scenario. A month after the burst, the flux is three orders of magnitude below that of 
the cosmological scenario, and, even if it were totally emitted in the soft X-ray band, totally unaccessible 
to current observational apparata. Thus an easy discrimination between the two models is possible. 

I now estimate the fraction of all flux radiated in the soft X-ray band, in the cosmological scenario. To 
do so, I have to resort to the customary treatment (Blandford and McKee 1977, Meszaros, Laguna and Rees 
1993, Sari, Narayan and Piran 1996) of electrons behind relativistic shocks. Thus I shall suppose that about 
half of the total internal energy is in a power-law distribution of electrons with index p = 2.5, giving rise to 
a synchrotron spectrum of index q = (p — l)/2 = 0.75 (Band et ai, 1993, Sari, Narayan and Piran 1996). 
Inverse Compton cooling does not contribute to emission in the relatively low energy bands I am interested 
in, so that the fraction fx of the total energy emitted in a band with photon energy ei < hv < Cu is 


fx = 


e 


(3-p)/2 _ i3-p)l2 

^W_ 

A3-p)/2 

fcm 


( 12 ) 


where Cm is the energy, in the observer’s frame, of the typical bremsstrahlung photons emitted by the 
highest energy electrons (ie., those having Lorenz factor 7 ^ in the shell’s frame): 




heB 

meC 


7m7 


(13) 


the Lorenz factor of the shell, 7 , appears here to convert the photons’ energy from the shell’s frame to the 
observer’s. I estimate 7 ^ by equating, as usual, the electron acceleration timescale with the synchrotron 
slow-down timescale, finding 7 ^ « 10^5“^/® (Meszaros, Laguna and Rees 1993). Inserting this into Eq. 

shows that Crn does not depend upon the magnetic field, and so it is independent of another major 
uncertainty of the problem, the efficiency with which equipartition magnetic fields are built up behind a 
relativistic shock. From Eq. and Eq. in the band 0.1 — 10 keV, for p = 2.5, I find 


fx 


0.2 

7I/4 


(14) 


Thus, the expected soft X-ray flux is 


Fx = 10 erg s ^ 


-2 l/3j^2/3 
cm Eg) 



(15) 
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where the function f(y)fx is also plotted in Fig. 1. The soft X-ray flux is of order « 10“^^ erg s~^ cm~^. 
The above equation shows why I concentrated on long bursts: short bursts having rj « 10^ (Sari and Piran 
1995) have soft X-ray fluxes which are lower than long bursts by about a factor of 20. This makes them 
unobservable with current satellites, even though still brighter than in the Galactic Halo model. 

A major uncertainty in estimating fx lies in postulating that the nonthermal electron population 
absorbs about half of all internal energy generated at the shock, but uncertainty in the slope of the spectrum 
also contributes. For p = 2, in fact, I would have found fx ~ O.l/y^/^, while, using = 10®, would have 
resulted in fx « O.Ol/y^^^. 

By analogy with supernova remnants, as the shock slows down, the fraction of total internal energy 
absorbed by the nonthermal electron distribution probably decreases, and becomes even harder to estimate. 
For this reason the theoretical prediction has not been extended to subrelativistic speeds of the shock. 


3. Discussion 

The most promising way to attack this problem is to try to identify afterglow emission from GRBs’ 
error box in the soft X-ray on a statistical basis. The advantage of doing this in the soft X-ray rather 
than in lower energy bands, where source contamination is also low, is that the expected fluxes from Eq. 

can be seen to be at least an order of magnitude higher. Higher-energy bands, instead, do not have the 
angular resolution necessary to avoid source confusion. In the radio band, where lower fluxes are offset by 
much larger collecting areas, a behaviour similar to that of the soft X-ray band is of course expected, but 
theoretical computation is made difficult by several subtleties. Consideration of this effect is thus postponed 
to a forthcoming paper. 

The small angular resolution of the WFCs onboard SAX (5 arcmin) allows follow up observations 
by narrow field instruments with the hope of little source confusion. This is exactly what has been done 
by Greiner et aL(1996), who identified three sources with the ROSAT HR Imager inside SAX error box. 
Of these, two sources are AGNs, while the third one, accounting for about half of the total flux (be., 

« 2 X 10“^® erg s~^ cm~‘^) detected by SAX and ASCA, has no optical counterpart, indicating a very 
unusual object. 

From Fig. 1, where the tickmark indicates the position of 43 post-burst days, and Eq. |^, we see that 
the expected flux, 5 x 10“^^ erg cm~^, compares remarkably well (and perhaps fortuitously, since we do 
not know the source distance) with the observation of 10“^® erg s~^ cm~^ (Piro et ai, 1996b, Murakami et 
al., 1996, Greiner et ai, 1996). Also, it should be noticed that Piro et aL(1996a) have set an upper limit to 
the soft X-ray flux from the burst region, immediately after the burst, of 10“^° erg s~^ cm~‘^ . Using Eq. 
and Fig. 1 it can be seen that the highest flux expected in this model is ^ 10“^^ erg cm~^. Thus, the 
current model reproduces correctly, for the most trivial choice of parameters, the observed features of the 
afterglow of GB960720. Hopefully, future observations ought to show that the afterglow has disappeared 
on a timescale of a few months, even though for this source it will be impossible to determine whether the 
light curve follows Fig. 1. 

Also, the apparent lack of interstellar absorption is consistent with the fireball model. The total baryon 
contamination « lO^^g (times at most a factor of 2 to include the mass swept up when the shell is still 
relativistic) when spread out over a spherical surface of radius Rgh < R < Rrei (Eqs. p| and ||) provides 
a total column depth Nh ~ 10^® — 10^^ cm~^, well within observational constraints. If furthermore. 
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as has been suggested, GRBs are related to mergers of neutron star binaries, they are expected to be 
distributed somewhat like pulsars, i.e. outside the galactic disk, where column depths do not approach the 
observational bounds. 

Lastly, I would like to point out the reason why detection of soft X-rays about a month after the 
burst can discriminate between cosmological and Galactic Halo models, since it is rather amusing. In 
the fireball model, there are 5 dimensional parameters: E, Mgj, pisM,c, D, where pisM is the density 
of the circumstellar matter, and D is source distance. Of these, c is a universal constant, and pisM 
is an external parameter, which ought to be considered as given, and which is, furthermore, relatively 
well-known, compared with the uncertainty in parameters such as E, Mej, each spanning several orders 
of magnitude. We can thus regard it as fixed. Thus, specifying that the two observational constraints, flux 
and time duration at Earth, be reproduced means fixing 2 of the 3 free parameters, leaving only one (say, 
D) undetermined. This corresponds to having a one-parameter (D) family of homologous solutions, each 
fitting observational data, each located at different distances from the observer. By making observations 
at a fixed time after the burst, about a month, we are observing cosmological and Galactic Halo models 
at non-homologous moments, thus breaking the similarity law that links them: in fact, the shock is still 
relativistic in the cosmological scenario, and well subrelativistic in the Galactic Halo scenario. 

In short, I have argued that, in the cosmological scenario of the hreball model, detectable fluxes (Eq. 
[isl and Fig. 1) of soft X-rays should be emitted in the two months following a gamma ray burst, with 
a non-thermal spectrum and a characteristic decrease (Fig. 1), while no detectable flux can arise in the 
Galactic Halo scenario. In particular, the expected flux of 10“^^ erg cm~‘^ compares remarkably well 
with the observations of GB960720 made about 40 days after the burst, by Piro et al.(1996b), Murakami 
et al., (1996) and Greiner et aL(1996). But the model also predicts the time-dependence of the afterglow 
(Fig. 1) and its disappearance after a few months, and is thus subject to more elaborate testing. Detection 
of the afterglow in a statistically meaningful sample of gamma ray bursts would simultaneously establish 
both the fireball model and the cosmological nature of GRBs. 


I am indebted to L. Stella and especially to Luigi Piro, for fruitful scientific conversations, and to an 
anonymous referee for constructive criticisms. 
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Fig. 1.— Solid line: plot of /(y) (from Eq. ^), against adimensional time ct^/Rsh where Rsh is defined in 
Eq. 1^, with r] = 100. Dashed line: plot of the product f{y)fx, where fx comes from Eq. |^. The tickmark 
represents, in scaled units, 43 days after the burst, the moment at which the follow up observations of GB 
960720 were carried out by Piro et al. (1996b). 
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